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Using Monte Carlo simulations, we study the character of the spin-glass (SG) state of a site- 
diluted dipolar Ising model. We consider systems of dipoles randomly placed on a fraction x of all 
sites of a simple cubic lattice that point up or down along a given crystalline axis. For x < 0.65 
these systems are known to exhibit an equilibrium spin-glass phase below a temperature Tsg oc x. 

At high dilution and very low temperatures, well deep in the SG phase, we find spiky distributions 
of the overlap parameter q that are strongly sample-dependent. We focus on spikes associated with 
large excitations. From cumulative distributions of q and a pair correlation function averaged over 
several thousands of samples we find that, for the system sizes studied, the average width of spikes, 
and the fraction of samples with spikes higher than a certain threshold does not vary appreciably 
with L. This is compared with the behavior found for the Sherrington-Kirkpatrick model. 

PACS numbers: 75.10.Nr, 75.10.Hk, 75.40.Cx, 75.50.Lk 


I. INTRODUCTION 

Complex systems are present in life and social sci¬ 
ences, information systems, and economics.® In these 
systems, different random distributions of their micro¬ 
scopic constituent parts give rise to diverse values of some 
macroscopic properties.121 A paradigmatic model in statis¬ 
tical physics that exhibits complexity is the Sherrington- 
Kirkpatrick (SK) model,® where the couplings between 
any pair of spins is randomly fixed to be ferromagnetic 
(FM) or antiferromagnetic (AF) regardless of the spin- 
spin distance. This model has both quenched spatial- 
disorder and frustrat ion ,® the two essential ingredients 
of spin glasses (SG).®^ Its exact solutioiP shows the 
existence of replica symmetry breaking (RSB)P differ¬ 
ent identical replicas of a sample JT with the same cou¬ 
plings may, in the thermodynamic limit, stay trapped 
in different states within the set of infinitely many pure 
states. These pure states are diverse, in the sense that 
they are sample-dependent. The distribution of the over¬ 
lap q between states of a given sample, Pj{q), is found 
to be a comb-like superposition of infinitely many 6- 
like spikes. In the macroscopic limit, after averaging 
over spatial disorder, the overlap distribution is given by 
p{q) = S{q — qra) + f{q) where f{q) is a non-zero smooth 
function for q < qm and zero otherwise. 

Whether the RSB picture describes correctly the be¬ 
havior of realistic spin glasses, such as dilute metallic 
alloys or concentrated insulator^ is still a matter of de¬ 
bate. The 3D Edwards-Anderson model (EA),li^in which 
only nearest-neighbor spins interact, is the simplest one 
with the essential ingredients of short-ranged SGs. No 
exact solution exists for the EA model, but there is con¬ 
sensus on the existence of a SG phase based on numerical 
simulations.!^ The applicability of a RSB scenario to the 
SG phase of the EA model is still controversial.!^^ In the 
so-called droplet picture, the SG phase is described in 
terms of a unique state (paired with the one obtained by 
a global spin inversion) with excitations that are com¬ 


pact droplets of the inverted state.!i^ According to this 
scenario, Pj(q) distributions do not exhibit diversity —in 
the sense that pj {q) is independent of J — and the aver¬ 
aged p{q) becomes a single delta function S{q — qm): p(q) 
is said to be trivial. Some trivial-non-trivial scenarios, 
between the droplet and RSB pictures, have also been 
proposed.!^! Early MG simulations point to a non-trivial 
scenario,!i^ but it has been found that the asymptotic 
behavior for p[q ~ 0) is only reached at very large sizes 
even for toy droplet models .!i^ 

There is growing interest in the study of sample-to- 
sample fluctuations of pj{q) from their average p{q)W^ 
Some recently proposed quantities give information on 
the heighfP^and average widtlP^of spikes found m.pj{q). 
By MG simulations, these quantities have been found 
to be nearly size-independent for the EA model, in con¬ 
tradiction with RSB predictions. However, these results 
have b een c riticized to be far away from the asymptotic 
regime.!2212il Some have found more useful to study the 
statistics of the area under pj{q) for q < The 

numerical study of SG models distinct from the SK and 
EA models has shed some light on the virtues and weak - 
nesses of these new probes for measuring diversity.l22l2ll 

Frustration in the SK and EA models comes from 
the competition between randomly distributed FM and 
AF couplings. However, frustration may also appear 
in fully occupied systems with no quenched disorder, 
such as the I sing model with pure AF interactions on 
a FCC lattice.!^ Dipoles packed in crystalline arrange¬ 
ments have frustration, exhibiting magnetic order that 
is strongly dependent on the lattice structure !^^!^^ ! Some 
ferroelectric^^ and magnetic insulators such as LiHoF 4 
are known to be well described by arrays of parallel Ising 
dipoles that behave as uniaxial ferromagnets.!^ In dipo¬ 
lar Ising systems (DIS), dilution put together with the 
built in geometric frustration results in SG behavior.!^ 
LiHoa;Yi_a;F 4 is one example that has been extensively 
studied. Experiment^ have found a SG phase for con¬ 
centration X = 0.16, and a FM phase for x > Xc where 
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Xc — 0.25. Recent MC simulations of systems of clas¬ 
sical Ising dipoles placed on a fraction x of the sites of 
the LiHoF 4 tetragonal lattice have found a SG phase for 
all 0 < cc < 0.25 for tempe ratures below a SG transition 
temperature Tsg oc 

Here we study a site-diluted system of dipoles, 
which are placed at random on a fraction x of the sites of 
a simple cubic (SC) lattice and point up or down along 
one of the principal axis. In the limit of low concentra¬ 
tions details of the lattice are expected to become irrel¬ 
evant. Therefore, our model at low concentrations has 
direct connection with the experimental and numerical 
work mentioned above. In previous MC worlP^ we have 
calculated the entire diagram of the system and found a 
SG phase for 0 < a: < ccc where x^ — 0.65 with the SG 
transition temperature given by Tsg{x) ~ x. We found 
from the following evidence that the SG phase behaves 
marginally: (i) the mean values qi =< |q| > decrease al¬ 
gebraically with L , (ii) averaged overlap distributions of 
q/qi appear wide and independent of L, (hi) ^l/T, where 
is a correlation length,l231 rises with L at constant T, 
but extrapolates to finite values as 1/L —?► 0. All of this is 
consistent with quasi-long-range order in the SG phase. 
Neither the droplet model nor a RSB scenario fit with 
this marginal behavior. 

The main aim of this paper is to study whether diver¬ 
sity may emerge in this geometrically frustrated model 
at low temperatures and high dilution, rather deep in its 
marginal SG phase by using the probes for measuring di¬ 
versity, in the sense that was specified above. The paper 
is organized as follows. In Sec. |TT]we define the model, 
give details on the parallel tempered Monte Carlo (TMC) 
algorithm,!^ and define the quantities we compute. We 
present results in Sec. m followed by concluding remarks 
in Sec. HV] 


II. MODEL, METHOD, AND MEASURED 
QUANTITIES 

A. Models 

We consider site-diluted systems of classical Ising spins 
on a SC lattice. All spins are parallel and point along the 
z axis of the lattice. At each lattice site a spin is placed 
with probability x. These spins are coupled solely by 
dipolar interactions. The Hamiltonian is given by 

H= T,ja,aj, ( 1 ) 

<*j> 

where the summation runs over all pairs of occupied sites 
i and j except i = j, Oi = ±1 on any occupied site i, 

T,j=ea{afr,gfil-izllrl), ( 2 ) 

where rij is the distance between i and j sites, Zij is the z 
component of r^, Sa is an energy, and a is the SC lattice 


constant. In the following, temperatures shall be given 
in terms of Sa/ks- 

Note that T^- values are not distributed at random, 
but depend only on the orientation of vectors Vij on a 
SC lattice. This i s why DIS exhibit AF order at con¬ 
centrations X > This is to be contrasted with 

random-axes dipolar models,l^in which Ising spins point 
along directions that are chosen at random, introducing 
randomness on bond strengths. 

In this paper we study DIS with x « Xc for which 
the details of the lattice structure are not important. Not 
surprisingly, in the limit of high dilution the behavior of 
DIS and the LiHoxY 4 _^F 4 system (a well known dipolar 
ferrom agnet for x = 1) have been found to be closely 

related. ^2133] 

For comparison, we study also the SK model: a set 
oi N = Ising spins Oi = ±I with interaction energies 
between any pair of spins at sites i and j given by JijUiGj 
with Jij = ±I/-\/]V chosen randomly, without bias, for 
all ij site pairs. 


B. Method 

For the models described in Sec. lH Al we have simulated 
a large number Ns of independent samples. By a sample, 
J, we mean a system with a given quenched distribution 
of empty sites for DIS (a quenched distribution of random 
couplings Jij for the SK model). The number of samples 
we average over is given in Tables I and H. We have 
tried not to make Ns smaller with increasing L. This is 
because statistical errors are independent of L, because 
of non-self-averaging. (However, for DIS with L = 10, 
we could only do 1.2 x 10"*^ samples. That took an Intel 
8-core Xeon processor E5-2670 some 60 years worth of 
GPU time). 

Thermal averages come from averaging over the time 
range [to,2to]i where to is the equilibration time. We 
further average over the Ns samples with different real¬ 
izations of quenched disorder. 

In order to accelerate equilibration at low tempera¬ 
tures in the glassy phase we use a parallel tempered 
Monte Carlo (TMC) algorithm.!^ We apply the TMC 
algorithm as follows. We run in parallel a set of n iden¬ 
tical replicas of each sample at different temperatures in 
the interval [Tmin,Tmax] with a separation A between 
neighboring temperatures. Each replica starts from a 
completely disordered spin configuration {ai}. We apply 
the TMC algorithm to any given sample in two stages. 
In the first stage, the n replicas of the sample N evolve 
independently for 10 MC Metropolis sweeps.l^ All dipo¬ 
lar fields throughont the system are ufated every time a 
spin flip is accepted. In the second stage, we give any 
pair of replicas evolving at temperatures T and T — A a 
chance to exchange states between them following stan¬ 
dard tempering rules which satisfy detailed balance.!^ 
These exchanges allow all replicas to diffuse back and 
forth from low to high temperatures and reduce equili- 
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bration times for the rough energy landscapes of SGs. 
We find it helpful to have the highest temperature Tmax 
larger than 1.6 x Tgg. We choose A such that at least 
30% of all attempted exchanges are accepted for all T. 

For DIS we use periodic boundary conditions (PBC). 
Details of the PBC scheme we use can be found in Ref.l^S. 
We let a spin on an occupied site i interact only with spins 
within an L X L X L cube centered on *. In spite of the 
long-range nature of the dipolar interaction, we do not 
perform Ewalds’s summations and exclude any contri¬ 
butions from repeated copies of the lattices beyond this 
box. This introduces an error which was shown for DIS 
in SC lattices to vanish as L —)■ oo, regardless of whether 
the system is in the paramagnetic, AF or SG phase (see 
Appendix I in Ref.l^. This result is not applicable to 
an inhomogeneous FM phase that may obtain on other 
lattices such as in LiHoF 4 . 


L 

T 

min 

T 

^ max 

A 

^0 

Ns 

4 

0.16 

1.60 

0.04 

lO'^ 

lO'^ 

6 

0.16 

1.60 

0.04 

lO'^ 

1.4 X 10^ 

8 

0.16 

1.60 

0.04 

2 X 10^ 

10’^ 


TABLE I. Simulation parameters for the SK model. The 
number of spins is N = L^, Tmin (Tmax) is the lowest (high¬ 
est) temperature and A is the temperature step in our TMC 
simulations. The number of MC sweeps for equilibration is 
to- Measurements are taken in the time interval [to, 2to]. The 
number of samples with different realizations of (quenched) 
disorder is Na. 


L 

T 

^ min 


T 

^ max 

Ai 

A 2 

^0 



Ns 

4 

0.05 

0.4 

1.65 

0.025 

0.05 

5 

X 

10® 

2 X 10® 

6 

0.05 

0.4 

1.65 

0.025 

0.05 

5 

X 

10 ® 

1.1x10® 

8 

0.05 

0.4 

1.65 

0.025 

0.05 

5 

X 

10 ® 

10® 

10 

0.075 

0.4 

1.65 

0.025 

0.05 

5 

X 

10" 

1.2x10^ 


TABLE 11. Same as in TableUbut for DIS with concentration 
X = 0.35. In our TMC runs for DIS we have chosen a tem¬ 
perature step of Ai for the temperature interval [TminjTm], 
and a bigger one, A 2 , for [Tm,Tmax]- 


C. Measured quantities 

Measurements were performed after two averagings: 
first over thermalized states of a given sample and second 
over a number of different samples. 

Given an observable u, we let wj = ( u )t stand for 
the thermal average of sample and u = [uj \j for the 
average over samples. 

We measure the Edwards-Anderson overlap 


parameter,!^ 

p) 

3 

where and are the spins on site j of identical 
replicas (1) and (2) of a given sample. Clearly, g is a 
measure of the spin configuration overlap between con¬ 
figurations of the two replicas. 

For each sample J we compute the overlap probability 
distribution pj{q). The mean overlap distribution p(q) 
over all replicas is defined by 

P{q) = [Pj{q)]j- ( 4 ) 

We also measure the mean square deviations of Pj{q), 
from the average p{q ), 

SpiQ^ = { {Pjil) - P{<l)V ]j- ( 5 ) 

In order to probe for RSB behavior we focus on over¬ 
laps between states that belong to different basins of 
atraction. With that aim, we compute the integrated 
probability functions defined by 

^ rQ 

^j= Pj{q)dq, ( 6 ) 

J-Q 

^J=if {Pj{<l)-p{Q)Vdqy^^, ( 7 ) 

J-Q 

and calculate their corresponding averages and . 
An advantage of working with quantities integrated over 
the interval q G {—Q, Q) is that statistical errors come 
smaller. 

Given that Xj is a (J^-dependent) random variable, it 
makes sense to explore how this variable is distributed. 
Following Reference |20], we define its cumulative distri¬ 
bution 11 ^( 2 :) as the fraction of samples having Xj < z. 

Yucesoy et have proposed very recently an ob¬ 
servable that is sensitive to spikes in the overlap distri¬ 
butions pj{q) of individual samples. They consider the 
maximum value of pj{q) for q G {—Q, Q), 

p^ = inax{ p^{q) :\q\ < Q }, (8) 

and count a sample as peaked if pj exceed some specified 
value. We compute the cumulative distribution nf(z) of 
Pj as the fraction of samples having pj < z. 

In previous paperJi^ we have obtained additional in¬ 
formation on the shape and width of spikes from a pair 
correlation function. Let fj{qi,q 2 ) = Pj{qi)pj{q 2 ), 

rQ rQ 

Gj{q) = / dqidq2 6{q2 - qi - q)fj{qi,q2), ( 9 ) 

^0 Jo 












4 


and let G^{q) be the average of Gj{q) over samples. We 
compute the normalized function 

g^{q)=G^{q)/ r G^{q)dq, (10) 

J-Q 

which is the conditional probability density that q = q 2 — 
qi, given that qi,q 2 & (0, Q). Note that g^{q) is largest at 
q^O, and that g^{q) = g^{-q), since pj(q) =pj{-q). 
It makes sense to define the width of g^{q) as 

vjQ = l/g^{0), (11) 

which is a measure of pattern thermal fluctuations for 
\q\ < Q. 

An additional interpretation of g^{q) is possible for 
sufficiently small T {T < OAT^g, roughly) so that indi¬ 
vidual spikes are clearly discernible. Assume, in addi¬ 
tion, that Q is sufficiently small so that contributions 
from samples with more than one spike in the 0 < q < Q 
range is negligible. Then, (i) finding on each sample one 
such spike, if there is one, (ii) calculating the self-overlap 
of such spike with a copy of itself shifted by a distance q 
(iii) adding the resulting function of q over all samples, 
and (iv) normalizing, gives g^{q). To that extent, g^{q) 
stands for an average over all spikes on the 0 < q < Q 
range. In Ref.l^^ we have also shown that if the width of 
spikes does not vary over different samples, then g^{q) 
is, for large systems, twice as wide as spikes are. 


D. Equilibration times 

We now explain how we make sure that thermal equi¬ 
librium is reached before we start taking measurements. 
To this end, some quantities are next defined. First, a 
pair of identical replicas of a given sample are allowed to 
evolve independently in time, starting at t = 0 from two 
uncorrelated random spin configurations. Let qt be the 
overlap between the configurations of the two identical 
replicas at time t. In addition, let 92 (t) be the average 
of q^ over all samples. During equilibration q 2 {t) is ex¬ 
pected to increase up to its equilibrium value. Semilog 
plots of q 2 {t) versus t displayed in Fig. for x = 0.35, 
L = 8,10 and the lowest used temperatures show that a 
stationary value <72 is reached only after some millions of 
MC sweeps. 

In order to check whether this stationary value q 2 is 
an equilibrium one, we define a second overlap, §t, not 
between configurations of pairs of identical replicas at the 
same time t, but between spin configurations of a single 
replica taken at two different times to and ti = to +1 oi 
the same MC run, 

qt{to) = N~'^'^aj{to)aj{to + t). (12) 

3 

Let q 2 {to,t) be the average of (qt(to))^ over all sam- 



FIG. 1. (Color online) Semilog plots of q2{to,t) and q 2 vs 
time t (in MC sweeps) for DIS systems with concentration 
X = 0.35 running at the lowest temperature for the values 
of L indicated in the hgure. T = 0.075 (T = 0.05) for L = 10 
(1/ = 8). q 2 is obtained from averages of q^ over time, starting 
at t = 0 from an initial random spin configuration. Here, 
to = 5 X 10^ MC sweeps for L = 10, and 5 x 10® MCS for 
L — 8. Data points at time t stand for an average over a time 
interval [t, 1.2t], and over 10® samples. 


pies. Suppose thermal equilibrium is reached long before 
time to has elapsed. Then, 52 (^ 0 , ^) and < 72 (t) should tend 
towards a common value <72 as t ^ to- Plots of q 2 {to,t) 
vs t are shown in Fig. for to = 5 x 10^ MC sweeps 
(5 X 10® MCS sweeps) for L = 10 (L = 8) for the same 
values of x and T as for (72 (t)- Note that both quantities, 
q 2 {to:t) and ( 72 (i); do become approximately equal when 
t —t to. In order to obtain equilibrium results, we have 
always chosen sufficiently large values of to to make sure 
that q 2 {to,t) « < 72 (t) for t >to- In our simulations, we let 
each system equilibrate for a time to and take averages 
over the time interval [to,2to]. All values of to and Ng 
are given in Table II. 

It has recently been shown that equilibration times in¬ 
crease with the roughness of the free-energy landscape 
of each individual sample.^ Numerous spikes in overlap 
distributions pj are the signature of samples that have 
numerous minima in their free-energy landscape. Visual 
inspection of overly distributions of samples like the 
ones shown in Fig. p shows fairly symmetric Pj curves 
even though some ofthem have several spikes. Then, our 
stringent equilibration criterion suggests that nearly all 
the samples are well equilibrated. 


III. RESULTS 

A. Overlap distributions 

As it has been found for other SG models, it is interest¬ 
ing to examine individual samples of DIS. In Fig.j^a) we 
plot pj(q) versus q for different samples at temperature 
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FIG. 2. (Color online) (a) Overlap distributions for DIS 

systems with L = 10, x = 0.35 and T = 0.1 for three samples 
with different realizations of disorder. For each sample we 
collect values of q over 5 x 10^ MC sweeps, (b) The same as 
in (a) but for T = 0.25 for the same sample set. Recall that 
the transition temperature is Tsg = 0.35. 



FIG. 3. (Color online) (a) Plots of the averaged distribution 
p[q) versus q for DIS systems with x = 0.35, T = 0.1, and 
the values of L shown, (b) Same as in (a), but for (Sp)^, the 
mean square deviations of pj{q)^ away from p{q) over all J 
samples. 


T/Tsg ~ 0.3. At this low temperature, some pj{q) dis¬ 
play well defined spikes centered on small q values, which 
seems to vary randomly from sample to sample. Qualita¬ 
tively similar di stributio ns have been observed for the EA 
and SK models.^ ^09 | 39 | jg clear that these inner peaks 
(that is, peaks away from g « ±1) come from overlaps 
between states that belong to different basins of attrac¬ 
tion. The main aim of this paper is to extract statistical 
information for these cross-overlap (CO) spikes situated 
on the interval q € (—QtQ)- Similar plots for higher 
temperatures (see Fig. I^b) for T/Tsg ~ 0.7), show that 
thermal fluctuations render individual spikes not clearly 
discernible. Then, in order to explore well within the 
SG phase, we have chosen the lowest temperature in our 
TMC simulations to be 0.2 Tgg. We report most of our 
results for T < 0.4 Tsg. We have also chosen a concentra¬ 
tion X = 0.35 which is far below the threshold for the AF 
phase. Both low temperatures and low concentrations 
result in large equilibrium times to- In addition, in or¬ 
der to obtain good sample statistics we need to simulate 
a large number Ns of samples for all system sizes stud¬ 
ied. All of this has restricted us to deal with relatively 
modest system sizes in our simulations. The simulation 
parameters are given in Tables I and II. 

Figure® a) shows the sample-averaged overlap distri¬ 
bution pi/^ for DIS at r = 0.1. At this low tempera¬ 
ture, p{q) exhibits two large peaks at ig™ with ~ 1 
and a relatively flat plateau with p(0) 0 in the region 

q S {—QtQ) for, say, Q « 1/2. The non-zero p{0) value 
does not change with L for the system sizes stud ied. This 
behavior, known for the SK and EA models,l^^is in con¬ 
tradiction with the droplet picture of SGs, for which p(0) 
vanishes as 


Plots of ((5p)^ vs p are shown for DIS at the same tem¬ 
perature in Fig. I^b). (Sp)'^, a measure of deviations 

of pj{q) from the average p{q), is clearly greater than 
p for g G (—1/2,1/2) indicating lack of self-averaging. 
More interestingly, (<5p)^ does not vary appreciably with 
L. This is at odds with the behavior found in MG 
simulations for the SK model for which (Sp)'^ oc L for 
T < 0.5 Tsg-^ Recall that in the RSB scenario, one ex¬ 
pects that Pj{q) exhibits with many sharp spikes in the 
region q G {—Q,Q) that become (5-hke functions as L 
increases, resulting in a diverging (dp)^ for macroscopic 
systems. 



FIG. 4. (Golor online) (a) Plot of versus T for DIS sys¬ 
tems with X = 0.35, Q = 1/2 and the values of L indicated 
in the figure, (b) The same plots for the SK model, with 
Q = 1/2 and the values of L shown. In both panels, error 
bars are smaller than symbol sizes. 
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B. Integrated overlap distributions 


Here we consider averages of both p and {5p)‘^ over 
q G {—Q,Q)- This allows us to focus on the contribu¬ 
tions of CO spikes and, in addition, to reduce statistical 
noise if Q is not too small. Plots of ^ the sample- 
averaged area under CO spikes, versus T are shown for 
Q = 1/2 in Figs, [djja) an dgb) for DIS and SK models re¬ 
spectively. In both cases, is, as far as we can see, size 
independent at temperatures well below Tgg. We obtain 
quantitatively similar results (not shown) for Q = 1/4. 
This is a strong piece of evidence against the validity of 
the droplet picture, for which X^ is expected to vanish 
as A similar behavior has b een fo und for the 

EA model in several MC simulationsHowever, it 
has been argued that strong finite-size effects may mask 
the asymptotic behavior at the system sizes currently 
available to MC simulation.l^^ Finally, we note that in 
Figs. I^a) and (b) X^ seems to vanish as T —>• 0 (as was 
long ago predicted for the SK model) .1^ 


Plots of vs T for DIS in Fig. j^a) show the pres¬ 
ence of finite size effects. Note that for small sizes, 
increases as T decreases only up to T = 0.15 (T = 0.075) 
for L = 4 (L = 6). We return to this point in Section 


HIC More interestingly, curves for larger sizes (L > 8) 
give a strong indication that does not diverge as L 
increases. This result is in contradiction with a RSB 
scenario and is in sharp contrast with the behavior ex¬ 
hibited in Fig. ^b) for the SK model, for which A*^ in¬ 
creases with yL at low temperatures. It is worth men¬ 
tioning that A*3 differs qualitatively from the average 
X 2 = [ {Xj)^ ]j when Q is not very small.^^ Ag has 
been investigated in detail in several papers and it is 
known to be size independent for both the EA and SK 
models.^ 



FIG. 5. (Color online) (a) Plot of Aq versus T for DIS systems 
with X = 0.35, Q = II2 and the values of L indicated in the 
figure, (b) Plot of versus T for the SK model, with 

g = 1/2. 



FIG. 6. (Color online): (a) Plots of vs q for the SK 
model at T = 0.2 with Q = 1/2 and the values of L shown, 
(b) Same plot as for (a) but for DIS systems with x = 0.35 at 
T = 0.075. (c) Same plot as for (b) but for T = 0.125. 


C. Pair correlation functions 


Data points for the pair correlation function for 
(5=1/2 are shown in Figj^a) and (b) for the SK model 
at T = 0.2 and DIS at T = 0.075 respectively. Note 
that these temperatures are such as T/Tsg ~ 0.2 in both 
cases. We find curves that are rather pointed with widths 
clearly smaller than Q. We obtain similar results for 
(5 = 1/4. Data for DIS in Fig. [^b) do not exhibit any 
significant size dependence. In contrast, g'^ curves for SK 
in Fig. [^a) become sharper as L increases. This result 
for the SK model is as expected for a RSB scenario. In 
the RSB solution, pj{q) is made of several cross-overlap 
spikes that for small values of q become 5-functions in 
the macroscopic limit, and densely fill the interval q € 
(—Q,Q). In striking contrast, our result for DIS suggest 
that the number of SG states do not grow with L at finite 
low temperatures. 

Some people have argued that comparing data for dif¬ 
ferent models (EA and SK) at the same value of T/Tsg 
is not meaningful. They find more appropriate to make 
such a comparison at temperatures for which X'^ values 
are the same.^ We follow this recipe and compare the 
data shown in Fig. [^c) for DIS at T = 0.125 and in 
Fig. [^a) for the SK model at T = 0.2. Apart from the 
fact that g^ becomes narrower as T decreases, we do not 
notice any qualitative difference. 

It is interesting to note that spikes, even though they 
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FIG. 7. (Color online) (a) Plots of wQ vs q for DIS systems 
with X = 0.35 for (3 = 1/2 and the values of L shown, (b) 
Plots of vs q for (a) but for the SK model, Q = 1/2 

and the values of L shown. In both panels all error bars are 
smaller than symbol sizes. 



z z 

FIG. 9. (Color online) (a) Plot of the cumulative distribution 
n{X®} versus x for DIS systems with x = 0.35, for T = 0.1 
and the values of L indicated in the figure, (b) The same 
plots as in (a) but for the SK model, for T = 0.16. 


have non-zero widths w at finite T, may not be discerned 
in pj{q) distributions of very small systems. Note that 
the minimum appreciable value of q for systems of N 
spins is given by Ag = 2/N. Thus, finite size effects are 
expected to come at very low T when w ^ Aq. This 
seems to be the case for the data shown in Figs.j^a) and 
[7];a) for L = 4 (L = 6) and T < 0.15 (T < 0.075^ 

Plots of versus T are shown in Figs. Ha) and Bb) 
for DIS and the SK model respectively, appear in 
Fig. [^a) to be size independent at least for L> 8. This 
points to finite widths for CO spikes in the L —?► oo limit 
for low (but finite) temperature. On the other hand 
values for the SK model displayed in Fig. [^b) appear to 
vanish as l/L^/"^ as L increases at least for T < 0.4, in 
agreement with the RSB picture. 



z z 

FIG. 8. (Color online) (a) Plot of the cumulative distribution 
n{p^} versus x for DIS systems with x = 0.35, for Q = 1/2, 
T = 0.1 and the values of L indicated in the hgure, where 
is the maximum value olpj{q) over the interval —Q < q < Q. 
(b) The same plot as in (a) but for the SK model, for Q = 1/2, 
T = 0.16 and the values of L shown. 


D. Cumulative distributions 

As interesting as they could be, pair correlation func¬ 
tions (PCFs) do only give information on how spiky sam¬ 
ple distributions Pj{q) are in the {—Q,Q) region. How¬ 
ever, PCFs do not give any information about the height 
of spikes located there. Following seminal work by Yuce- 
soy et fl/.li^we study here UP{z), the fraction of samples 
without any spike in {—Q,Q) with height larger than z. 
Plots of UP versus z for DIS at T = 0.1 are shown in 
Fig. I^a). They give a strong indication that reach 
a size-independent shape for L > 8 for a wide range of 
values of z. In order to check for the robustness of our 
UP{z) values, we have grouped all available samples in 
K ensembles of 10^ samples each, calculated UP for each 
ensemble k = 1,K, and obtained the standard devia¬ 
tion (SD) of the K resulting values. Tiny vertical bars 
in Figs. 1^ and stand for such SDs. Plots in Fig. ^a) 
are to be compared with the ones displayed in Fig. ^b) 
for the SK model. Note in Fig. [^b) that, at least for 
z > 0.5, UP(z) clearly decreases as L increases, indicat¬ 
ing a proliferation of high spikes for larger sizes. This 
is as expected for the RSB picture, for which CO spikes 
become (5-like functions in the thermodynamic limit. 

Finally, we report results on how the random variable 
is distributed. Previous MC simulations have found 
a similar behavior for the EA and SK models when deal¬ 
ing with the cumulative distribution n^(z) of quantity 
z = Xj The mean-field theory of the SK model of¬ 
fers precise predictions on H^ (z) and their moments.^ 
In particular, n^(z) is found to follow a power law for 
small z. For small values of U^{z) oc z^, where 
y stands here for . Some people have found it use¬ 
ful to study the median of the cumulative distribution, 
which is predicted to reach a non-zero value in the ther¬ 
modynamic limit in the RSB picture but vanishes for 
the droplet model. Log-log plots of U^ versus X for 














DIS displayed in Fig. [^a) show curves with power-law 
behavior for small z. Data do not show any significant 
deviation for sizes L > 8. The median (marked by the 
crossings points of the curves with the horizontal dot¬ 
ted line in the figure) decreases as L increases reaching a 
non-zero value. The counterpart plots for the SK model 
are shown in Fig. [^b). In agreement with previous MC 
work on the SK modelpSEH we find strong finite-size ef¬ 
fects. However, curves seem to converge to some limiting 
curve as L increases.!^ Note that, in contrast with DIS, 
the median of H^ increases as L increases. All our re¬ 
sults for for both the SK model and DIS are not in 
contradiction with a RSB scenario. 


IV. CONCLUSIONS AND DISCUSSION 

By tempered Monte Carlo calculations, we have stud¬ 
ied the low temperature behavior of a diluted system of 
classical dipoles placed on a SC lattice. These dipoles are 
Ising spins randomly placed on a fraction x of all lattice 
sites and point up or down along a common crystalline 
axis. 

Previous MC studieJi^ for this model have provided 
strong evidence for the existence of a SG phase for 
X < 0.65 with a SG transition temperature Tsg{x) ~ x. 
The SG phase was then found to have quasi-long range 
order, as in the 2D-XY model.l^ Neither the droplet 
model nor a RSB scenario fit with this marginal behav¬ 
ior. Despite the existence of this soft SG order, we find 
in our simulations spiky overlap distributions pj {q) that 
are strongly sample-dependent, as previo usly fou nd in 
simulations for the EA and the SK models.^^^ttHall 

We have studied the statistics ofpj{q) for q G {—Q, Q) 


using some recently proposed observables 
that p{q) and Sp{q) (as well as their integrated counter¬ 
parts and A'3) do not vary appreciably with L. 

From a suitable defined pair correlation functiorP^ we 
compute an averaged width that appears to remain 
finite as L increases. Complementary to this result, we 
find that the fraction of samples with spikes higher than 
a certain threshold does not vary appreciably with L. All 
of this points to finite width for CO spikes in the L —>■ oo 
limit at low temperatures. 

Our results are in clear contradiction with droplet 
model predictions. On the other hand, a direct compari¬ 
son of our data with MC data obtained for the SK model 
shows that crucial RSB predictions are also at odds with 
some of our results. It is noteworthy that the findings 
enumerated above for DIS are strikingly similar t o the 
ones found in previous MC work for the EA model,^ ^^° l ^^ l 
though finite size effects have been reported to be strong 
for the EA model for the available systems sizes in MC 
simulations. 
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